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The inclusive production of ui and <f> mesons is studied in the backward region of the interaction of 
12 GeV protons with polyethylene, carbon, and copper targets. The mesons are measured in e + e~ 
decay channels. The production cross sections of the mesons are presented as functions of rapidity 
y and transverse momentum pr- The nuclear mass number dependences (A dependences) are found 

tO be J 40' 710 ±0-021(stat)±0.037(syst) for u mes(ms and A 0.937±0.049(Btat)±0.018( ByB t) for ^ megons fa ^ 

region of 0.9 < y < 1.7 and p T < 0.75 GeV/c. 
PACS numbers: 13.85.Ni, 13.85.Qk, 14.40.Cs 



I. INTRODUCTION 

The modification of the vector meson spectral function 
in hot and/or dense matter is currently a hot subject in 
terms of spontaneous chiral symmetry breaking and par- 
tial restoration of the symmetry in nuclear matter. Cur- 
rently, this is the focus of several experiments Q, 0, 0, Q ■ 
The experiment KEK-PS E325 was performed to mea- 
sure the vector meson spectral functions in dense matter, 
i.e., nucleus. Thus far, we have reported the signature of 
the mass modification of vector mesons 0, IS 13 ■ These 
observations are fairly remarkable; hence, we also per- 
formed analyses to determine the absolute cross sections 
and nuclear mass number dependences of the production 
of these mesons in order to understand the underlying 
production mechanism. 

The nuclear mass number dependence of the cross sec- 
tions for the particle production is usually parameterized 
as 



a(A) = a A a 



(1) 
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for a target nucleus with mass number A. When the 
collision energy ^/snn is sufficiently large, the parameter 
a in the production of light mesons like pions or p mesons 
is about 2/3 HQ. 

This can be interpreted by considering 
such productions to be dominated by primary collisions 
at the front surface of a target nucleus. Note that the 
mean free path of incident protons in a nucleus is as small 
as 1.4 fm. On the other hand, a tends to unity in the 
case of hard reactions like the production of J/tp at high 
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energies, ^snn > 20 GeV 

The present experiment was performed at \J s^vw = 
5.1 GeV. At higher energy, a for <fi meson production was 
reported to be 0.81 ± 0.06,0.96 ± 0.04, and 0.86 ± 0.02 at 
T^vw = 11.6 E3, 14.2 and 15.1 GeV Q, respec- 
tively. However, there is no reason to believe that these 
values are applicable at our energy. A few heavy ion in- 
duced experiments at ^/snn from 4.9 to 5.4 GeV 0,0] 
reported <f> meson production data, to which the present 
experiment can provide complementary data. Note that 
there have been no measurements of w mesons and <p 
mesons with p + A reactions at snn around 5 GeV. 

A (j) meson is almost a pure ss state. Therefore, 
the production of a meson without other accompany- 
ing strange particles is suppressed by the Okubo-Zweig- 
Iizuka rule [13. This results in the effective threshold 
energy being as high as ^/snn = 3.9 GeV, which corre- 
sponds to 2m p + 2m ss -, since two ss pairs are effectively 
needed to realize the <f> production. Our collision energy 
is fairly close to this effective threshold. In addition, 
the importance of additional mechanisms for <f> meson 
production, such as an intrinsic ss component in nucle- 
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FIG. 1: Schematic view of the experimental setup from 
the top, which is designed symmetrically with respect to the 
beam. The components of this setup are referred as left and 
right arms in this article. 

oris and 4>pn coupling 0, |2(j , have been suggested 
by theorists; however, they have thus far been insuffi- 
ciently studied at our energy from an experimental view- 
point. Thus, to understand the production mechanism at 
our energy, basic measurements such as production cross 
sections and a parameters are indispensable. 

In this article, we present inclusive production cross 
sections and a parameters of u> and <p mesons measured 
via the e + e~ decay channels in 12 GeV p+p, p + C, and p 
+ Cu collisions. The results are compared to those from 
nuclear cascade simulation JAM [21| , and the implications 
to production mechanisms are discussed. 



II. EXPERIMENT 

The spectrometer was built in the EP1-B beam line 
at the 12 GeV proton synchrotron in KEK. The beam 
line was designed to deliver a 12 GeV primary proton 
beam with an intensity of up to 4 x 10 9 protons per 
spill. The beam was extracted for 1.8 s with a repeti- 
tion rate of 1/4 Hz. The beam intensity was monitored 
with about 10% accuracy using ionization chambers [2^ 
located downstream of the spectrometer. 

Figures ^and [3 show a top view and a side view of the 
experimental setup, respectively. The spectrometer was 
designed to simultaneously measure e + e~ and K + K~ 
pairs. A dipole magnet and tracking devices were com- 
monly used together with electron and kaon identification 
counters. 

The magnet was operated at 0.71 T at the center of 
the dipole gap and provided 0.81 Tm of field integral 
for the tracking. The magnetic field map was calculated 
by the finite-element-analysis software TOSCA j^. The 
calculated map agreed well with the measured map, and 
the difference was negligible when compared to the mo- 



FIG. 2: Cross section of the experimental setup along the 
plane of the center of the kaon identification counters. 



mentum resolution of this spectrometer. During the data 
collection periods, the field strength was monitored every 
4 s with a nuclear-magnetic-resonance (NMR) probe lo- 
cated at the center of the surface of the lower pole piece. 
The magnetic field map was scaled run by run according 
to the NMR data. The fluctuation of the magnetic field 
was found to be less than 10~ 5 within a typical run of 
two hours long. 

The targets were aligned inline along the beam axis 
at the center of the magnet. The target materials and 
configurations are shown in Table [I] 

The beam profile in the horizontal direction was mea- 
sured by counting the interaction rate by changing the 
beam position at the target. The beam position was 
moved by the bending magnet located about 10 m up- 
stream of the target. In this measurement, the center 
target, whose thickness was 1 mm, was rotated by 90° 
around the vertical axis, and was used as a 1 mm- wide 
probe. The typical beam size in the horizontal direction 
was found to be 2.0 mm in full-width half-maximum. The 
beam size in the vertical direction was known to be al- 
most the same as in the horizontal direction as seen in 
a view of a luminescence plate which was temporarily 
inserted during the beam tuning. 

Three tracking devices — the vertex tracking chamber 
(VTC), cylindrical drift chamber (CDC), and barrel drift 
chamber (BDC) — were used to determine the trajectories 
of the charged particles. In the present analysis, the mo- 
mentum was determined using CDC and BDC. The mo- 
mentum resolution a p was a p = (1.37% • p) 2 + 0.41% 2 • 
p (GeV/c), where p is the momentum of a particle. 

For the electron identification, two types of gas 
Cerenkov counters (FGC and RGC) and three types of 
lead glass calorimeters (SLG, RLG, and FLG) were em- 
ployed. The gas Cerenkov counters were horizontally 
segmented into 6°. The radiator of the gas Cerenkov 
counters was isobutane at room temperature and atmo- 
spheric pressure. The refractive index is 1.00127 at the 
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TABLE I: Summary table for the targets, beams, and trigger modes used in the present analyses. In the 2002 run, four Cu 
targets were used. 



year 


target 


position 
(mm) 


interaction 
length (%) 


radiation 
length (%) 


number 
of protons 


trigger 




CH 2 


-48 


0.111 


0.195 




electron 


1999 


C 


38 


0.106 


0.213 


3.0 xlO 13 


/ kaon 




Cu 


-7 


0.0391 


0.412 






2002 


C 





0.213 


0.431 


3.2 xlO 14 


electron 




Cu 


±24, ±48 


4 x 0.0539 


4 x 0.565 







standard temperature and pressure, which corresponds 
to a momentum threshold of 2.26 GeV/c for charged pi- 
ons. All lead glass detectors were built with SF6W |2J| 
and were typically segmented to 3.5° horizontally. The 
typical energy resolution for 1 GeV electrons is 0.15 GeV. 
The acceptance for electrons ranged from ±12° to ±90° 
in the horizontal direction and from —22° to 22° in the 
vertical direction. 

The counters for the kaon identification were start tim- 
ing counters (STC), hodoscopes (HC), aerogel Cerenkov 
counters (AC), and time of flight counters (FTOF). The 
acceptance for kaons ranged from ±12° to ±54° in the 
horizontal direction, and from —6° to 6° in the vertical 
directions. The three counters (STC, AC, and FTOF) 
were horizontally segmented to 6°, and HC was typically 
segmented horizontally to 3°. The time of flight (TOF) 
of charged particles was measured using STC and FTOF 
with a resolution of 0.36 ns and flight length of about 
3.7 m. Kaons and pions were separated in a momentum 
range from 0.53 to 1.88 GeV/c using an aerogel with a 
refractive index of 1.034. 

The electron trigger signal required a hit in FGC ac- 
companied by a geometrical coincidence with RGC, SLG, 
or RLG. To select events containing e + e~ pairs with large 
opening angles, both the left and the right arms were re- 
quired to contain at least one e + or e~ candidate. The 
typical efficiency of the trigger for the electron pairs in 
the acceptance was 92.4%. 

The kaon trigger signal was obtained from the coinci- 
dence of STC, HC, and FTOF. The charged pion contam- 
ination was reduced by using AC as a veto trigger. Pro- 
ton contamination was reduced by setting a TOF window 
for the kaons by using a rough momentum value calcu- 
lated by combining the hits in STC, HC, and FTOF in 
the trigger. 



III. DATA ANALYSIS 

A. Analysis of lo and 4> — > e + e~ in p + C and p + 
Cu collisions 

1. Event reconstruction 

The charged tracks were reconstructed from hit posi- 
tions in the drift chambers by using the Runge-Kutta fit- 
ting method. After the tracks were reconstructed, tracks 
corresponding to momenta between 0.4 and 2.0 GeV/c 
were selected for further analyses. The lower limit in the 
momentum range was set based on the threshold of the 
trigger, whereas the upper was set in order to avoid pion 
contamination above the gas Cerenkov threshold. 

Pairs of positive and negative tracks were required to 
satisfy the trigger condition. All the e + and e~ candi- 
dates were reexamined so that an FGC hit association 
could be obtained with an RGC, SLG, or RLG associa- 
tion, depending on the location of the track. For candi- 
dates associated with lead glass calorimeters, the momen- 
tum ratio E/p should be larger than 0.5 for energy in the 
calorimeters to be obtained. We chose the value of 0.5 
balancing the purity of e + and e~ with the statistics of 
the present data. Figure [3] shows the distribution of the 
energy and momentum of the present data with FGC as- 
sociations. In this figure, it is clearly seen that electrons 
range along the line of E = p. After an e + e _ pair was 
detected, we simultaneously refit the e + and e~ tracks 
by constraining them to have the same vertex point on 
the interaction target. Finally, we identified 5.69 x 10 5 
e + e~ pairs. 1 

The reconstruction efficiency of the tracks from the 
targets was evaluated by both an eye-scan and a detector 
simulation using GEANT4 j2(J . 

In the former method, first we used only one of the 
two arms and determined the target in which an interac- 
tion occurred. Then we visually scanned about 200 event 
displays and found track candidates in the other arm us- 



The number of recorded events for electron data were 
7.41 x 10 7 and 5.08 x 10 8 in 1999 and 2002, respectively. 
A detailed description of the spectrometer can be found 
in Ref. I^. 



In a further analysis, we used only events in which e+ went into 
the left arm and e~ went into the right arm because the number 
of events satisfying the opposite criterion comprised only 6% of 
the data. 
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FIG. 3: A contour plot of the distribution of energy and mo- 
mentum of the present data with FGC associations. The solid 
line shows the criterion of electron identification described in 
this article. 

ing the drift chamber hit information with help from the 
interaction point. The tracking efficiency was evaluated 
by seeing if the tracking program could find those eye- 
scanned track candidates. The efficiency was found to be 
67% on average. 

In the detector simulation method, the reconstruction 
efficiency was evaluated using simulation tracks embed- 
ded onto the real data. According to this method, the 
efficiency was found to be 78%. 

In order to combine the results of the two evaluations, 
we simply assumed the average 73% to be the efficiency 
in the present analysis and considered half the difference 
as a systematic error. The reason of this discrepancy 
is unknown, and this is one of the major sources of the 
systematic uncertainty in the present data analysis. 

The inefficiency of the vertex reconstruction process 
was evaluated as follows. We took all the combinations of 
e + and e~ tracks regardless of the position of the closest 
point of each pair, and obtained the yield of the u> meson 
peak. Then we used the vertex reconstruction program 
and the required the event vertex belong to any of the 
target disks. In addition, we fit the e + and e~ tracks 
together by constraining them to have the same vertex 
point on the interaction target, and required that the x 2 
over the number of degrees of freedom (NDF) should be 
less than 5. As a result, we lost 6.0% of the ui meson yield, 
so that the vertex reconstruction efficiency was 94.0%. 

Contaminations due to the misidentification of pions 
and other particles in the present data were evaluated 
at the mass region of the u> meson, i.e., from 0.75 to 
0.80 GeV/c 2 , by tightening the electron identifications 
with gas Cerenkov counters and lead glass calorimeters 
until the misidentification becomes negligibly small. In 



this mass region, in the p + C data, we found that 18% 
of the events result from misidentification and 18% are 
from uncorrelated e + e~ pairs. 



2. Corrections 

Besides the tracking efficiency, several detector effects 
were evaluated and corrected as described below. The ef- 
ficiencies of the electron identification counters were eval- 
uated as a function of momentum by using pure electron 
samples from 7 conversions and Dalitz decays. These 
electron samples were identified as a zero-mass peak in 
the e + e~ spectra, and they were not required to partic- 
ipate in the trigger to avoid trigger bias. The obtained 
efficiencies were typically 85% for FGC, 86% for RGC, 
and 97% for the lead glass counters. 

The energy losses of the tracks through the detec- 
tors were estimated by using GEANT4 simulation. Typi- 
cally, the reconstructed momentum gave a value lower by 
3 MeV/c for a 1 GeV/c electron due to the energy loss. 
The momentum difference was corrected track by track 
by using a correction table obtained by the simulation. 
It should be noted that this correction compensated only 
mean energy loss. For effects that cause an eventual large 
energy loss like bremsstrahlung, the corrections were car- 
ried out in a different manner, as described later. 

The geometrical acceptances for vector mesons V were 
obtained as functions of the invariant mass, rapidity y, 
and transverse momentum px by the simulation. The 
acceptances were averaged over azimuth ip, and isotropic 
decays of V — > e+e - were assumed. The effects of the 
trigger, i.e., requirements of a geometrical correlation be- 
tween the electron identification counters, were also con- 
sidered. In order to obtain the yields of vector mesons, 
mass spectra were corrected for the acceptance in the 
mass range above 0.55 GeV/c 2 . Below 0.55 GeV/c 2 , the 
correction was too large to evaluate reliable values. The 
obtained acceptances at the o> and the 4> meson masses 
are tabulated in Fig. 0] 



3. Spectrum decomposition 

In Fig.|SJ the left and the right panels show the spectra 
of the invariant mass of e + e _ pairs in the range of 0.9 < 
y < 1.7 and px < 0.75 GeV/c without and with the 
acceptance correction, respectively. The peaks of to and 
4> are distinctly observed. 

The e + e~ mass spectra were fit and decomposed into 
the dielectron decays u> — > e + e~, <j) — ► e + e~, and 
p° — > e + e~; the Dalitz decays 77 — ► 7e + e _ , u> — > 7r°e + e~, 
<p — ► 7r°e + e~, and <j> —* rje + e^\ and the combinato- 
rial background. The origins of the combinatorial back- 
ground were pairs which were picked up from two inde- 
pendent Dalitz decays or 7 conversions, and pairs like 
e^Tr^ due to the misidentification. The Dalitz decay 
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FIG. 4: Acceptances (%) for lo meson mass (left) and <f> meson mass (right) as functions of y and pr of the mesons. The errors 
are statistical and are obtained from the simulation. 



7T° — > 7e + e~ contribution is negligible in the acceptance 
of the present data. 

For the invariant mass distributions of the lo — > e + e~ 
and <f> ~~ > e + e~ decays, we used in the fit a Breit-Wigner 
function 

XX = H (2) 

dm ( m - mo ) 2 + r t 2 ot /4 

convoluted with a Gaussian function for the experimental 
resolution. If the relativistic Breit-Wigner shape is used 
instead, the results do not change significantly. Here, a 
is the cross section; m, the invariant mass of the e + e~ 
pair; N, a normalization factor; mo, the meson mass; 
and r to t, the natural decay width. For the shape of the 
p° — > e + e~ decay, we used the relativistic Breit-Wigner 
shape 



dm ( 9 2 \ 2 , 9-n2 

\rn z — m p0 \ +m^l tot 

instead of the Breit-Wigner function. To obtain the mass 
distribution in the spectrometer acceptance, we used the 
GEANT4 simulation with an input momentum distribu- 
tion of p° mesons that was obtained using JAM. The 
e + e _ invariant mass spectra from the Dalitz decays of 77, 
lo, and cj) mesons were also obtained by the simulation. 
Their e + e~ distributions were determined by following 
the vector meson dominance model given in Ref. |27| by 
using the mother meson distributions obtained by jam. 
The combinatorial background shape was evaluated us- 
ing an event mixing method, combining e + and e~ tracks 
picked from different events. 

The free parameters of the fit were the yields, the 
peak positions, and the mass resolutions of lo and cf> 
mesons; the yields of r\ and p° mesons; and the num- 
ber of the background events. As mentioned earlier, 
the spectra were corrected in a mass range only above 
0.55 GeV/c 2 . Therefore, below 0.55 GeV/c 2 , the uncor- 
rected spectra were used in the fit mainly to accurately 



obtain the amount of the background. Although the fit 
region was from 0.20 GeV/c 2 to 1.2 GeV/c 2 , the mass 
regions from 0.600 to 0.765 GeV/c 2 and from 0.955 to 
0.985 GeV/c 2 were excluded in order to avoid the effect 
of the excesses below the lo and cf> peaks, which in other 
publications |H S we have claimed as the signal of 
the mass modification. In the present analysis, however, 
we did not assume any underlying physics for the excess, 
and we aimed only to obtain the yields of lo and <p mesons 
correctly, as discussed later. 

The fit results are shown by thin solid lines in Fig. [S] 
The x 2 / NDF were obtained as 157/152 and 192/152 for 
the p + C and p + Cu interactions, respectively. 

After obtaining the raw yields of the lo and </> mesons 
by the fit, the correction for hard energy losses such as 
bremsstrahlung or a radiative tail was applied. The hard 
energy loss causes a low mass tail in the invariant mass 
distribution. The loss of the yield due to this tail could 
not be evaluated using the procedure described above; 
hence, we performed studies using the GEANT4 simula- 
tion. The yields of the tails with respect to the integrals 
of the Breit-Wigner peaks were found to be 11.9 ± 1.0% 
for the lo mesons and 10.4±0.8% for the </> mesons. These 
values were simply added to the peak yields. The uncer- 
tainty of these corrections includes an ambiguity with 
respect to the cross sections of such hard energy losses in 

GEANT4. 

By the fit procedure described above, the peak posi- 
tion and mass resolution of the <p meson were obtained as 
1019±1 MeV/c 2 and 11.8X1.0 MeV/c 2 , respectively. The 
peak position of the <f> meson is consistent with the values 
of the Particle Data Group [2^] , and the mass resolution 
is consistent with the simulation value of 10.7 MeV/c 2 ; 
however, the peak position of the lo meson was found to 
be lower by 2.2 MeV/c 2 . The width of the omega mesons 
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FIG. 5: The e + e~ mass distributions of the acceptance-uncorrected data (left) and the corrected data (right) in the kinematic 
range of 0.9 < y < 1.7 and pr < 0.75 GeV/c. The lines in the left panels represent the backgrounds, Dalitz components, p° 
components, and their sum. The amount of e + e~ decays of to and 4> mesons is represented by thick solid lines in the right 
panels, and the thin solid lines represent the sum of all the components. The number of <j> meson Dalitz decays is negligibly 
small. 



was also broader than that expected by the simulation. 2 
Here, we adopted a conservative approach in order to 
obtain the yield of unmodified lo mesons by estimating 
the additional systematic error due to this peak shift. 
The error for the u> meson yield was evaluated as 0.48% 
for a p + C interaction and 1.8% for a p + Cu interac- 
tion, which were obtained by forcing the peak position 
in the fit to the higher value of 2.2 MeV/c 2 . The mo- 
mentum scale and resolution were also verified using the 
K s — ► 7r + 7r~ data 0. The measured peak position of K s 
was 496.8 ± 0.3 MeV/c 2 , which was consistent with that 
of the simulation result, 496.9 MeV/c 2 . The measured 
mass resolution of K s was 3.9 ± 0.4 MeV/c 2 , which was 
expected to be 3.5 MeV/c 2 by the simulation. 



2 If we take internal bremsstrahlung into account, the peak po- 
sitions of both the u> mesons and the <f> mesons agree with the 
expected positions. 



4- Systematic errors 

In addition to the uncertainties described previously, 
the following systematic uncertainties were studied to ob- 
tain the cross sections. 

The uncertainty in the background shapes could be 
a source of the systematic error. The background was 
a result of uncorrelated pairs that were obtained from 
two independent Dalitz decays, 7 conversions, or other 
meson decays. The background shape was obtained by 
an event mixing technique. In the mixing process, it is 
possible that the correlated e + e~ pairs from the decays 
of p°, ui, and 4> mesons deform the estimated background 
shape; this could result in a systematic error. In order to 
estimate the systematic error, two methods were used in 
the event mixing. One was to use all electrons except for 
those belonging to the oj and (f> meson mass regions, which 
are 0.765 to 0.800 GeV/c 2 and 0.995 to 1.035 GeV/c 2 , 
respectively. In the other method, we used all the pairs 
in the event mixing but with weights in order to obtain 
a self-consistent shape of the background. The weights 
were obtained as a function of the e + e~ mass as the 
ratio of the background shape to the real mass spectrum, 
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and they were self-consistently determined by repeating 
the fit several times. Wc adopted the latter method to 
determine the background shape. The difference between 
the two methods was assigned as a systematic error. The 
difference in the yield of to mesons was 0.23% and that 
in <fi mesons was 0.60%. 

The normalization and spectral shape of the combina- 
torial background were affected by all the other corre- 
lated e + e~ components since they were obtained by the 
fit. The systematic errors in the peak-yield determina- 
tion due to the ambiguity of the Dalitz decays of rj, to, 
and <f) mesons were evaluated by doubling or eliminating 
those yields and refitting the mass spectra. These sys- 
tematic errors were found to be 0.62% and 0.85% for u> 
and 4> mesons, respectively. 

It was difficult to consider the uncertainty of the p° 
shape, since the mass region from 0.600 to 0.765 GeV/c 2 
could not be represented by known sources, and the shape 
of the mass modification was not well understood. In 
this analysis, we simply considered the relativistic Breit- 
Wigner shape for the p° distribution and performed a fit 
by excluding the excess region. In order to evaluate the 
systematic error, we fixed the p° yield at zero as an ex- 
treme case and reperformed the fit. In this fit, the u> yield 
increased by 3.19% for the carbon target and 7.65% for 
the copper target; further, the <f> meson yield increased 
by 3.45% and 5.48% for each of the above targets, respec- 
tively. We considered these values as systematic errors 
due to the unknown p° distribution. It should be noted 
that the u) mesons can also be modified such that the to 
cross sections obtained in the present analysis are only 
for that component whose shape is consistent with the 
un-modified shape. 



TABLE II: Systematic uncertainties for ui and 4> meson yields 
for the data taken in 2002. 





beam intensity 




10 % 




C target thickness 




0.28% 




Cu target thickness 




0.55% 




trigger electronics 




3.8 % 




track reconstruction 




14.9 % 




vertex reconstruction 




4.7 % 




acceptance 


1.23% 




1.25% 


hard energy loss 


1.18% 




0.94% 


electron identification 


2.92% 




2.60% 


background shape 


0.23% 




0.60% 


mass scale for p + C 


0.48% 






mass scale for p + Cu 


1.83% 






Dalitz yield 


0.62% 




0.85% 


p° from p + C 


3.2 % 




3.5 % 


p° from p + Cu 


7.6 % 




5.5 % 



total for p + C 


19.5 


% 


19.5 


% 


total for p + Cu 


20.8 


% 


20.0 


% 



TABLE III: Systematic uncertainties for a parameters of u> 
and d> mesons for the data taken in 2002. 



acceptance 


1.47% 


1.13% 


target thickness 


0.52% 


0.39% 


background shape 


0.27% 


0.54% 


mass scale 


1.13% 




Dalitz yield 


0.37% 


0.48% 


p° yield 


4.89% 


1.25% 


total 


5.2 % 


1.9 % 



A possible cause of another error might lie in the effi- 
ciency estimations for electron identification. These un- 
certainties were evaluated by using the error bands of 
efficiency curves of electron identification counters and 
were obtained as 2.92% for cu mesons and 2.60% for <fi 
mesons. 

All the systematic errors are summarized in Table iHl 
In the evaluation of the absolute cross section, the errors 
of the beam intensity, target thicknesses, and efficiency 
of the trigger electronics were also considered. In sum- 
mary, the systematic errors for the cross sections for to 
mesons were 19.5% for the carbon target, and 20.8% for 
the copper target, and those of <f> mesons were 19.5% and 
20.0%, respectively. 

The systematic errors for the a parameters are only 
due to the items that differ between carbon and copper- 
targets. These are listed in Table ITTT1 In summary, the 
uncertainties for the a parameters of to and <fi mesons 
were 5.2% and 1.9%, respectively. 



B. Extraction of the cross section in p + p 
interaction 

The target subtraction method was employed to ob- 
tain the cross section in p+p collisions; this was achieved 
by using the data taken in 1999. The production cross 
sections <j(j>) in p + p collisions were calculated by the 
formula a(p) = (ct(CH2) — er(C))/2 in the region of 
0.9 < y < 1.7 and pt < 0.6 GeV/c. The analysis was 
performed in almost the same manner as that for the 
data in 2002. 

To minimize the uncertainty due to experimental dif- 
ferences between the 1999 and 2002 data, we normalized 
the yield of to mesons with carbon targets in the 1999 
data to that in 2002 data and extracted the absolute 
cross sections in p + p collisions. Another uncertainty 
arose from the difference in the target thicknesses that 
were measured with an accuracy of 0.8%. We obtained 
20.1% and 20.0% as the systematic uncertainty in p + p 
collisions for the lo meson production and <j) meson pro- 
duction, respectively. 
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TABLE IV: Branching fractions in the tables of the Particle 
Data Group H|. 



decay mode 



branching fraction 



e e 
K+K 



(6.96 ±0.15) x 10" 
(2.96 ±0.04) x 10" 



49.2 



+0.6 



% 



IV. RESULT AND DISCUSSION 



A. Global feature of the obtained results 

From the decay branching fractions listed in Ta- 
ble IIVI the differential cross sections of the inclu- 
sive u> meson production are obtained as 14.30 ± 
0.34(stat)±2.79(syst) mb-c/GeV and 46.63 ± 1.25 (stat)± 
9.70(syst) mb-c/GeV in p + C and p + Cu colli- 
sions in the region of 0.9 < y < 1.7 and pr < 
0.75 GeV/c, respectively; further, those of <f> meson pro- 
duction are 0.270 ± 0.017(stat) ± 0.053(syst) mb-c/GeV 
and 1.290 ±0.070(stat)±0.258(syst) mb-c/GeV, respec- 
tively. 3 The a parameters of u> and <j) mesons are ob- 
tained as 0.710 ± 0.021 (stat) ± 0.037(syst) and 0.937 ± 
0.049(stat) ± 0.018(syst), respectively. The difference is 
0.227 ± 0.054(stat) ± 0.041(syst), and it is statistically 
significant. 

In order to compare the production cross sections in 
the e + e~ decay channels with those in the K + K~ de- 
cay channel and the previous measurement of 12 GeV/c 
p + p — > p°X by Blobel et al. 29], the present data 
were extrapolated to the backward hemisphere — region 
of xf < or y < 1.66 — where xf is Feynman's x. The 
production cross sections in the backward hemisphere are 
listed in Table [V] The correction factors for this extrap- 
olation from the measured regions — 0.9 < y < 1.7 and 
Pt < 0.75 GeV/c — to the backward hemisphere were cal- 
culated by using JAM, since the shapes of pt and y spec- 
tra are consistent with the result of JAM calculation, as 
described later. 

The cross sections of the inclusive 4> meson production 
in the K + K~ decay channel are also listed in Table [V] 
These were obtained from the previous analysis in this 
experiment [30l l3l| . The a parameter of the 4> meson 
production measured in the K + K~ decay channel is ob- 
tained as 1.01 ± 0.09 using the p + p, p + C, p + Cu, and 
p + Pb data in the spectrometer acceptance; this value 
is statistically consistent with the present e + e~ analysis. 

Figure[S]shows the cross sections in the backward hemi- 
sphere as a function of the nuclear mass number. The 
dotted lines represent the result of the JAM calculation 




1 10 10 2 

nuclear mass number 

FIG. 6: Production cross section in the backward hemisphere 
as a function of mass number. The circles, open squares, and 
filled squares show the cross section of U) mesons, </> mesons 
measured in the e + e~ decay channel, and those in the K + K~ 
decay channel, respectively. The vertical lines represent the 
statistical errors, and the brackets represent the systematic 
errors. The previous p + p — > p°X measurement |29j is indi- 
cated by a triangle. The dotted lines represent the prediction 
by JAM simulation results. The solid and dashed lines repre- 
sent the q parameterization (see text). 



and the solid and dashed lines represent the measured a 
parameterization. The a parameters shown in the figure 
were 0.710, 0.937, and 1.01 for the data of u> — > e + e~, 
<t> — ► e + e~, and <j) — > K + K~, respectively. It should 
be noted that these values were obtained from the data 
within the spectrometer acceptance before they were ex- 
trapolated to the backward hemisphere. 

The previous measurement |2^| yielded the total cross 
section of 1.8 ± 0.25 mb for p° mesons. By assuming a 
p° /w ratio of 1.0 ± 0.2, 4 we obtain the u> meson produc- 
tion in the backward hemisphere as 0.90 ± 0.22 mb for a 
comparison with the present measurement. The triangle 
in Fig. represents the obtained value. The measured cj 
meson production cross section in the present p + p colli- 
sion data is consistent with their p° cross section within 
the error. 



3 In this article, neither the effect of internal bremsstrahlung 
the uncertainty of the branching fractions is considered. 



4 The p°/u> ratio was measured by the reactions of p 4 
charged particles and p + p — > u> + charged particles 
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TABLE V: Meson production cross section in the backward hemisphere. The first errors are statistic, where as the second 
errors are systematic. The data points are plotted in Fig. |5] 





uj (mb) 


(e+e") (mb) 


(K + K~) (mb) 


H 


2.3 ± 1.1 ± 0.5 


0.034 ± 0.045 ± 0.007 


0.03 ± 0.19 ± 0.01 


C 


13.4 ± 0.3 ± 2.6 


0.240 ± 0.015 ± 0.047 


0.39 ± 0.05 ± 0.19 


Cu 


49.0 ± 1.3 ± 10.2 


1.21 ± 0.07 ± 0.24 


1.84 ± 0.27 ± 0.89 


Pb 






6.0 ± 1.8 ± 2.9 



B. Differential cross section of uj and 4> production 
measured in e + e~ decays 

In a region 0.9 < y < 1.7 and p T < 0.75 GeV/c, the 
differential cross sections of uj and (j) mesons were ob- 
tained for each y or px bin, as shown in Fig. These 
are also listed in Table IVT1 

The previous p + p — » p°X measurement [29j was plot- 
ted using triangles, as shown in Fig. It should be 
noted that the data points were scaled by factors of 5.81 
(= 12 71 ) for clarity. 

The distributions obtained by the JAM calculation are 
compared with the measurements in Fig. [7] The to- 
tal cross sections obtained from the JAM calculation are 
larger than that obtained from the data. Hence, they 
are scaled by the factors 0.489 for p + C — > ujX, 0.421 
for p + Cu -> ujX, 1.006 for p + C -> <f>X, and 0.686 
for p + Cu — > 4>X. These scale factors were determined 
as the ratio of the total cross sections in the acceptance 
between the data and JAM. The shapes of the differen- 
tial cross sections of JAM are consistent with the present 
data, although the absolute cross sections are systemat- 
ically larger than the data. 

Figure [H] shows the a parameters of u> and <f> mesons 
obtained in the region of 0.9 < y < 1.7 and px < 
0.75 GeV/c. These are also listed in Table EH The 
flat lines represent the averaged a parameters with the 
errors shown in the left column. The dotted and dashed 
curves represent the results of JAM for p + C and p + Cu 
collisions. Although the a parameters of the JAM calcu- 
lation are significantly larger than those of the data, the 
difference between uj and <f> mesons in the a parameters 
is similar to what is seen in the JAM calculation. 



C. Discussion 

An interesting characteristic is the difference in the y 
dependence between uj and <f> meson production cross sec- 
tions. While the uj production is essentially independent 
of y, the <fi production increases towards the smaller y val- 
ues. In addition, the a parameters of uj and <f> meson pro- 
duction are different by 0.227± 0.054(stat) ± 0.041(syst). 
This significant difference in the a parameters confirms 
that the production mechanisms of uj and <j) mesons are 
different. 

One possible mechanism of <f> meson production is a 



hard reaction between incident and projectile nucleons. 
If the hard reaction is dominant like J/ip production at 
higher energies, a is expected to be around unity and 
independent of y. Further, the y distribution of the cross 
section is expected to be symmetric with respect to y c . m ., 
which is not the case in the present <j) production data. 

Another possible effect to explain the observed char- 
acteristics of the a parameter and cross sections is the 
effect of secondary collisions in a target nucleus. In this 
case, no hard reaction is necessary. These effects are ex- 
pected to increase in a smaller rapidity region for p + A 
interaction, and the cross sections and a are expected to 
be larger in the backward region. 

Although the JAM does not reproduce the data quan- 
titatively, it yields fairly similar shapes in the inclusive 
cross sections of uj and </> meson productions; further, it 
also predicts the difference between the a parameters of 
these mesons. In the jam calculation, more than 90% of 
the 4> mesons are produced by secondary collisions mostly 
between non-strange mesons and nucleons in a target nu- 
cleus, while uj mesons are produced both in primary and 
secondary reactions. For example, 30% and 50% of the uj 
mesons are produced in secondary reactions in the cases 
of 12 GcV p + C and p + Cu reactions in jam, respec- 
tively. These differences in the production mechanism in 
JAM will surely be an important step to understand the 
measured differential cross sections. 

The difference between the scaling factors applied to 
</> meson production in p + C and p + Cu collisions, in 
jam, described in Sec. II V Bj can be understood qualita- 
tively by the overestimated contribution of the secondary 
collisions in jam. Although in those calculations almost 
all the 4> mesons are produced in secondary collisions, the 
present data suggest that the contribution of the primary 
collisions can be larger. For the case of the uj meson, jam 
just predicts larger absolute production yields than seen 
in the present data. 

In summary, we measured the inclusive differential 
cross sections of the uj and <fi meson production in p + 
A collisions in the backward region. The difference in 
the a parameters between uj and (j) mesons confirms that 
the production mechanisms of uj and (f> mesons are dif- 
ferent. The results are compared to the nuclear cascade 
calculations. 
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FIG. 7: Differential cross sections of uj (top) and (f> (bottom) mesons as functions of y (left) or pt (right). The statistical 
errors are represented by vertical bars, and the systematic errors are represented by brackets. The previous p + p — > p°X 
measurement |2S| is indicated by triangles after scaling up with the nuclear mass number dependence obtained in the present 
analysis. The dashed and solid curves represent scaled cross sections of JAM in p + C and p + Cu collisions, respectively. The 
rapidity of the center of mass system is 1.66, as indicated by y c .m.- The data points are listed in Tabic IvTI 
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: o data, co 
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0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 0.2 0.4 0.6 0.8 1 1.2 

y y c.m. p T (GeV/c) 



FIG. 8: The a parameters of u) (circle) and <f> (square) mesons as functions of y (left) or px (right). The statistical errors are 
represented by vertical bars and the systematic errors are represented by brackets. The horizontal lines represent the averaged 
a parameters of u> and <j> mesons with the errors indicated on the left side. The dotted and dashed curves represent the a 
parameters of 4> an d w mesons by JAM. The rapidity of the center of mass system is indicated by y c .m.- The data points are 
listed in Table EH 
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TABLE VI: Differential production cross sections in 0.9 < y < 1.7 and pr < 0.75 GeV/c measured via e + e decay channel. 
The first errors are statistical and the second are systematic. The data points are plotted in Figs. Q and |H1 

uj from C uj from Cu 



(mb c/GeV) (mb c/GeV) 





0.9 -1.0 


14.7 


± 1.5 


± 2.9 


47.8 


± 4.6 


± 10.0 


0.708 ± 0.083 ± 0.042 




1.0 -1.1 


14.2 


± 1.2 


± 2.8 


51.2 


± 4.2 


± 10.7 


0.770 ± 0.070 ± 0.041 




1.1 -1.2 


15.0 


± 1.1 


± 3.0 


48.8 


± 4.6 


± 10.2 


0.706 ± 0.071 ± 0.040 


y 


1.2 -1.3 


14.5 


± 0.8 


± 2.8 


48.3 


± 2.9 


± 10.1 


0.724 ± 0.049 ± 0.039 




1.3 -1.4 


14.4 


± 0.8 


± 2.8 


45.8 


± 2.7 


± 


9.6 


0.693 ± 0.048 ± 0.039 




1.4 -1.5 


14.4 


± 0.7 


± 2.8 


45.3 


± 3.2 


± 


9.5 


0.689 ± 0.052 ± 0.039 




1.5 -1.6 


15.3 


± 0.7 


± 3.0 


46.9 


± 2.9 


± 


9.8 


0.671 ± 0.045 ± 0.038 




1.6 -1.7 


14.2 


± 0.7 


± 2.8 


46.7 


± 3.8 


± 


9.7 


0.716 ± 0.056 ± 0.039 




0.00-0.12 


5.7 


± 0.4 


± 1.1 


18.8 


± 1.2 


± 


3.9 


0.714 ± 0.057 ± 0.038 




0.12-0.24 


15.1 


± 0.7 


± 3.0 


49.2 


± 3.2 


± 10.2 


0.709 ± 0.047 ± 0.038 


PT 


0.24-0.36 


18.1 


± 0.8 


± 3.5 


63.1 


± 3.4 


± 13.2 


0.750 ± 0.042 ± 0.039 


(GeV/c) 


0.36-0.48 


19.5 


± 0.9 


± 3.8 


64.2 


± 3.6 


± 13.4 


0.714 ± 0.044 ± 0.038 




0.48-0.60 


16.8 


± 1.3 


± 3.3 


52.6 


± 3.6 


± 11.0 


0.687 ± 0.063 ± 0.040 




0.60-0.75 


13.9 


± 1.1 


± 2.7 


41.3 


± 3.9 


± 


8.6 


0.655 ± 0.075 ± 0.041 


0.9 < y < 1.7, p T < 0.75 GeV/c 


14.30 


± 0.34 


± 2.79 


46.63 


± 1.25 


± 


9.70 


0.710 ± 0.021 ± 0.037 








<f> from 


C 




4> from C 


u 




a 




0.9 -1.1 


0.348 ± 0.046 ± 0.068 


1.60 


± 0.17 


± 


0.32 


0.916 ± 0.101 ± 0.022 


V 


1.1 -1.3 


0.232 ± 0.032 ± 0.045 


1.33 


± 0.13 


± 


0.27 


1.050 ± 0.101 ± 0.020 




1.3 -1.5 


0.277 ± 0.029 ± 0.054 


1.20 


± 0.11 


± 


0.24 


0.881 ± 0.084 ± 0.020 




1.5 -1.7 


0.255 ± 0.037 ± 0.050 


0.93 


± 0.13 


± 


0.19 


0.780 ± 0.119 ± 0.019 


p T 


0.00-0.25 


0.185 ± 0.024 ± 0.036 


0.93 


± 0.10 


± 


0.19 


0.971 ± 0.101 ± 0.019 


(GeV/c) 


0.25-0.50 


0.405 ± 0.032 ± 0.079 


1.78 


± 0.13 


± 


0.36 


0.890 ± 0.066 ± 0.019 




0.50-0.75 


0.255 ± 0.032 ± 0.050 


1.19 


± 0.16 


± 


0.24 


0.924 ± 0.111 ± 0.021 


0.9 < y < 1.7, pr < 0.75 GeV/c 


0.270 ± 0.017 ± 0.053 


1.290 ± 0.070 ± 


0.258 


0.937 ± 0.049 ± 0.018 
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